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Locating  Unimoduiar  Roots 


Introduction 

The  well-known  formula  of  de  Moivre  ([4],  pg.  22),  implies  the  polynomial  q(z)  = 
z7  —  1  has  seven  roots  symmetrically  distributed  on  the  unit  circle  with  one  at  z  =  1 
and  the  others  equally  spaced  at  angular  intervals  of  2tt/7.  In  contrast,  inserting  a 
lower  order  term  as  in  p(z)  =  z7  +  z5  -  1  results  in  a  less  symmetric  pattern  of  the 
roots  as  shown  in  Figures  1  and  2.  Observe,  however,  that  two  of  the  seven  roots  of  p(z) 
still  appear  to  lie  on  the  unit  circle.  This  is  indeed  the  case  as  we  shall  see  below. 


Figure  1.  The  unit  circle  and  roots  of  Figure  2,  The  unit  circle  and  roots  of 

q(z)=z7-  1.  p{z)  =  z7  +  £-1. 

Changing  the  middle  term  of  p(z)  from  z5  to  z4  results  in  the  polynomial  having  no 
roots  on  the  unit  circle.  Furthermore,  if  we  consider  the  family  of  polynomial  equations 
z7  +  zk  —  1  =  0  for  k  =  1, 2, , , . ,  6,  then  the  only  case  in  which  there  are  solutions 
of  unit  modulus  is  when  k  =  5.  This  observation,  along  with  other  similar  examples, 
prompted  the  investigations  which  led  to  this  paper. 

We  will  refer  to  roots  that  lie  on  the  unit  circle  as  unimoduiar  roots.  We  ask  when 
the  polynomial  p (z)  —  zTt  +  zk  ~  1  has  unimoduiar  roots,  and  provide  an  answer  in 
terms  of  conditions  on  n  and  k ,  Here  we  assume  that  n  and  k  are  integers  with  n  >  2 
and  1  <  k  <  n  -  1.  Throughout  this  paper,  we  let  C  denote  the  unit  circle  in  the 
complex  plane. 

Inserting  the  lower  order  term  zk  in  the  polynomial  q(z)  =  zn  —  1  results  in  a 
complicated  pattern  for  the  location  of  the  roots.  For  example,  Figure  3  shows  all  380 
roots  for  the  family  of  polynomials  given  by  p(z)  —  z20  +  zk  —  1}  where  1  <  k  < 
19,  It  appears  that  some  of  these  roots  are  unimoduiar,  but  most  are  not.  We  will  see 
that  when  k  —  4,  eight  of  the  twenty  solutions  to  z20  +  z4  —  1  =  0  are  unimoduiar, 
but  for  k  7^  4,  there  are  no  unimoduiar  solutions. 

For  a  different  example,  if  n  =  17,  then  the  corresponding  family  of  polynomials, 
p{z)  =  z1 7  F  zk  —  1.  has  exactly  two  unimoduiar  roots  when  k  —  1 }  7,  or  13.  For  all 
other  values  of  k  there  are  none. 

Studying  the  roots  of  trinomials  is  an  old  problem  dating  back  to  the  nineteenth 
century.  Mel  man's  recent  paper  [2]  provides  improved  bounds  on  both  angular  sectors 
and  annular  regions  containing  the  roots.  It  also  gives  historical  references  to  several 
of  the  original  papers  on  this  topic. 
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Figure  3>  The  unit  circle  and  the  380  roots  of  p(z)  =  z20  4  zk  -  1,  for  !  <  k  <  19. 


More  generally,  the  problem  of  counting  the  number  roots  of  a  polynomial  that 
occur  inside,  on,  or  outside  the  unit  circle  has  been  extensively  studied.  The  Schur- 
Cohn  procedure  (see  [3])  provides  a  recursive  method  for  counting  these  quantities 
for  a  given  polynomial  However,  the  trinomials  we  study  here  are  singular  cases  for 
the  Schur-Cohn  algorithm.  There  are  variations  of  the  algorithm  to  handle  the  singular 
cases,  but  they  are  difficult  to  implement  and  unlikely  to  yield  an  explicit  formula 
for  the  number  of  uni  modular  roots  in  terms  of  n  and  k.  Another  shortcoming  of  the 
Schur-Cohn  approach  is  that  the  algorithm  fails  to  specify  the  exact  locations  of  the 
roots. 

In  our  analysis,  we  are  able  to  characterize  the  number  and  locations  of  uni  modular 
roots  of  p(z)  =  zTt  4-  zh  —  I  in  terms  of  conditions  on  n  and  k. 


Necessary  and  Sufficient  Conditions 

Let  p(z)  =  zn  4  zk  —  1  and  assume  p  has  a  uni  modular  root  e%e .  Substituting  elB  into 
the  equation  p(z)  =  0  and  rearranging  gives 

+  =  L  (1) 

The  only  pair  of  unimodular  complex  numbers  that  sum  to  1  are  1/2  ±  iyf  3/2-  This 
follows  by  observing  the  pair  must  be  conjugates  and  are  constrained  to  satisfy  the 
equation  of  the  unit  circle.  In  polar  form,  the  numbers  1  /2  ±  i x/3/2  can  be  written  as 

Focusing  our  attention  on  the  angles  nQ  and  kQ>  we  see  there  exists  integers  a  and 
/?  satisfying  the  pair  of  equations 


nQ  =  4-  27ra 

o 

(2) 

k$  =  +  2*fi. 

(3) 
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Solving  both  equations  for  8  and  equating  the  results  gives 

„  it  27TQ  7T  27 t3 

±  4.  _ -  =  q; - (_  - 

3  n  n  3  k  k 

Multiplying  through  by  3nk/ir  and  simplifying  leads  to 

n  +  k  —  ±6(n/3  —  ka).  (4) 

The  preceding  discussion  establishes  the  following  lemma: 

Lemma  1.  If  p(z)  =  zn  +  zk  -  1  has  unimodular  roots,  then  n  4-  k  =  0  mod  6. 

The  result  in  Lemma  1  gives  a  necessary  condition  on  n  and  k,  but  the  condition 
is  not  sufficient.  Indeed,  the  polynomials  z4  +  z2  -  1  and  z°  +  z3  -  1  both  satisfy 
n  +  k  =  0  mod  6,  but  neither  has  any  unimodular  roots.  We  observe  in  these  two  ex¬ 
amples  that  gcd(n,  k)  ^  1.  Making  the  substiutions  w  =  z'2  and  u  =  z3  lead  to  the 
associated  polyomials  w2  +  w  -  I  and  u3  +  u  -  1  respectively.  Neither  of  these  sim¬ 
pler  polynomials  have  any  unimodular  roots.  This  observation  motivates  the  following 
lemma: 

Lemma  2.  Let  n,  k  €  N,  with  1  <  k  <  n  -  1,  and  let  g  =  gcdfn,  k).  Then  p(z)  = 
zn  q -  zk  _  l  has  unimodular  roots  if  and  only  if  g(z)  =  zn/s  +  zh/a  -  1  has  unimod- 
ular  roots. 

Proof.  Assume  A  is  unimodular  and  thatp(A)  =  0.  Then  Xs  is  also  unimodular  and 

?( =  (X3)n/s  +  (X3)k/a  -  1  =  Au  4-  Afc  -  1  =  p( A)  =  0. 

Conversely,  assume  7  is  unimodular  and  3(7)  =  0.  Let  u  represent  any  of  the  uni¬ 
modular  gih  roots  of  7;  so  w9  =  7.  Then 

p(u)  =  (w)n  +  (w)fc  -  1  -  (u9)n/s  +  (u9)k/<>  -  1  =  9(7)  =  0. 


Lemma  2  makes  it  clear  that  we  only  need  to  consider  polynomials  p{z)  =  zn  + 
zk  -  1  in  which  n  and  k  are  relatively  prime.  In  cases  where  gcd(n,  k)  =  g  >  1,  we 
instead  consider  the  polynomial  q(z)  =  zn/,J  +  zk/g  -  1.  The  unimodular  roots  of  p 
are  in  g-to-1  correspondence  with  the  unimodular  roots  of  q. 

With  the  restriction  gcd(n,  k.)  =  1,  we  can  state  and  prove  the  converse  to 
Lemma  1.  The  proof  will  depend  on  the  following  classical  result  from  the  theory 
of  linear  Diophantine  equations  (see  [1],  pg.  35  for  a  complete  discussion). 

Theorem  1.  The  integer  equation  ax  +  by  -  c  has  a  solution  if  and  only  if  g  divides 
c,  where  g  =  gcd(ct,  b).  Moreover,  if  x  =  Xq  and  y  =  y0  is  one  particular  solution, 
then  the  complete  integer  solution  is  given  by 


,  0  a 

x  =  x<j  +  rn-  y  =  y0-m-j 

9  9 

where  rn  G  Z. 

Lemma  3.  If  n  -f-  k  =  0  mod  6  and  gcd(n,  k)  =  1,  then  are  the  only  uni¬ 
modular  roots  of  p{z)  =  zn  +  zk  —  1. 
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Proof  The  hypotheses  on  n  and  k  imply  that  neither  is  divisible  by  2  or  3  (exercise). 
Thus  both  integers  are  of  the  form  ±1  mod  6.  Hence,  there  exist  non-negative  integers 
$  and  t  such  that  n  =  6s  ±  I  and  A;  =  6t  =F  1- 

We  assume  n  =  6 s  +  1  and  k  =  6f  -  1,  with  the  other  case  being  very  similar. 
Thus,  we  can  write 

p(z)  =  z6*+1  +  z*"1  -  1  -  z“*z  +  Z6t  Q)  -  1. 

A  direct  calculation  shows  p  (e±£^/3)  =  0. 

Now  assume  p{ei&)  —  0,  We  will  show  that  8  =  ±|  up  to  addition  by  a  multiple 
of  2tt>  We  note  that  equation  (4)  was  derived  under  the  assumption  that  p  had  a  uni- 
modular  root.  Equation  (4)  is  really  two  linear  diophantine  equations  in  the  variables 

a  and  /J : 


— Qka  +■  6  nfi  =  n  +  k  (5) 

6 ka  -  6 n/3  —  n  4-  k.  (6) 

To  see  that  Theorem  1  applies  we  need  gcd(6A:,  6 n)  to  divide  n  +  fc.  It  readily  follows 
from  gcd(n7  k)  —  1  that  gcd(6&.  6rc)  —  6,  and  we  assumed  n  4-  k  is  divisible  by  6. 

Recalling  that  n  =  6s  +  1  and  h  —  6t  -  1,  we  leave  it  as  an  exercise  for  the  reader 
to  verify  that  the  pair  a  =  s,  /3  =  t  is  a  solution  to  equation  (5),  and  the  pair  a  — 
n  -  s,  jS  =  k  -  t  is  a  solution  to  equation  (6),  Thus,  by  Theorem  1,  the  complete  set 
of  integer  solutions  to  equation  (5)  is  given  by 

a  =  s  +  m  «  n  /?  =  t  +  m  *  ky  (7) 

where  m  €  Z.  Similarly,  the  complete  set  of  integer  solutions  to  equation  (6)  is  given 
by 

a  =  n  -  s  -  m  ■  n  ft  =  k  -  t  -  m  A\  (8) 


where  m  G  Z. 

Now,  adding  equations  (2)  and  (3),  making  the  substitutions  72.  —  Gs  4- 1  and  k  — 
-  1,  and  solving  for  8  yields 


77  a  +  /3 
3  s  t 


(9) 


Substituting  the  possible  values  for  a  and  /?  from  (7)  and  (8)  into  equation  (9)  gives 
8  =  mod  2tt  as  required,  ■ 


The  Complete  Result 

We  can  now  completely  characterize  when  and  where  the  imimodular  roots  occur  for 
our  family  of  polynomials. 

Theorem  2.  Let  p(z)  =  zn  +  zk  -  1  and  let  g  —  gcd(n,  k ).  If  6  divides  n/g  4 
k/g7  then  p  has  exactly  2 g  imimodular  roots ,  zm  and  zmy  given  by 


",  (  7T 

27TtoY 

=  exp 

i  1%  + 

;  JJ 
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where  0  <  m,  <  g  —  1 . 

As  an  illustration  of  the  theorem,  suppose  that  n  =  70.  By  trial  and  error  we  see 
that  there  are  only  six  values  of  k  satisfying 

70  k 

gcd(70,  k)  +  gcd(70.  *)  ”  °  m°d 
These  values  and  the  number  of  unimodular  roots  are  shown  in  Table  1. 


k 

gcd(70,  k) 

Number  of  unimodular  roots 

2 

2 

4 

14 

14 

28 

26 

2 

4 

38 

2 

4 

50 

10 

20 

62 

2 

4 

Table  L  For  n  —  70,  the  values  of  k  resulting  in  unimodular  roots. 


In  the  preceding  example,  there  were  six  A-values  for  which  the  resulting  polyno¬ 
mial  had  unimodular  roots.  There  are  64  unimodular  roots  out  of  4830  total  roots  for 
the  entire  family  p(z)  =  z70  +  zk  -  1.  The  number  of  unimodular  roots  for  a  family 
varies  widely  with  the  value  of  n.  In  constrast  to  the  case  where  n  =  70,  there  are 
arbitrarily  large  values  of  n  such  that  for  all  1  <  k  <  n  -  1  there  are  no  unimodular 
roots;  any  n  —  2**  ■  3*%  where  ^  v  €  N  is  an  example. 

The  case  where  n  >  5  is  prime  is  more  tractable  .  Here,  the  hypotheses  of  Theorem  2 
are  satisfied  by  the  smallest  integer  kQ  for  which  n  +  k0  is  a  multiple  of  6.  Starting 
from  k0y  there  are  also  unimodular  roots  for  h  =  kQ  +  Cm,  m  =  0, 1 , , . . ,  as  long  as 
k  <  n  —  1 .  In  each  such  case*  there  will  be  exactly  one  pair  of  conjugate  unimodular 
roots. 

Finally*  the  requirement  that  n  4-  k  =  0  mod  6,  with  gcd(n.  A;)  —  1  implies 
that  n  —  5  and  k  —  1  are  the  smallest  possible  values  of  n  and  A;  that  result  in 
unimodular  roots. 


Problems  for  Investigation 

We  conclude  this  paper  by  describing  possible  areas  for  further  research.  Some  of  these 
questions  may  be  suitable  for  advanced  undergraduate  projects. 

1 .  For  fixed  n,  what  is  the  total  number  of  unimodular  roots  for  the  family  of  polynomi¬ 
als  p(^)  =  zn  +  zl  —  1,  where  1  <  k  <  n  —  1?  It  is  straightforward  to  answer  this 
question  when  n  is  prime.  However,  as  shown  in  the  example  above  with  n  =  70 T  it  is 
more  difficult  when  n  is  composite. 

2.  An  obvious  extension  of  our  problem  is  to  insert  two  lower  order  terms.  Thus,  consider 
the  family  of  polynomials  given  by  p(z)  =  zn  +  zk  +  zj  -  1,  where  1  <  j  <  k  < 
n  —  1  Special  cases  such  as  j  —  k  -  1  and  j  =  n  —  k  might  good  places  to  start. 

3.  For  all  pairs  n  and  A’,  we  observed  that  roots  occur  both  inside  and  outside  the  unit 
circle.  It  is  a  straightforward  application  of  Rouche's  theorem  (see  [4]}  to  determine 
an  annular  band  about  C  in  which  all  the  roots  must  lie,  and  to  see  that  the  band 
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gets  narrower  with  increasing  values  of  n.  Fixing  n  and  allowing  k  to  vary,  numerical 
experiments  show  complicated  patterns  for  the  number  of  roots  inside  or  outside  [he 
unit  circle.  However,  if  one  views  half  of  the  unimodular  roots  as  occuring  inside  and 
half  outside  the  unit  circle,  then  more  predictable  patterns  occur.  Restricting  to  the  case 
where  n  is  prime  and  letting  k  increase,  we  find  the  number  of  roots  inside  or  outside 
the  unit  circle  as  a  function  of  k  only  changes  when  n  H-  k  =  0  mod  6, 

4*  As  seen  in  Figure  3t  the  set  of  roots  from  the  family  of  polynomials  p(z)  —  zn  - b 
zk  —  1  (with  n  fixed)  presents  a  striking  image*  It  is  interesting  to  animate  this  process 
by  plotting  the  complete  set  of  roots  for  each  value  of  k^  as  k  increases  from  1  to 
n  —  L  The  animation  shows  that  the  roots  appear  to  loop  around  the  nth  roots  of 
unity  at  different  rates*  We  are  happy  to  provide  our  Mathematica  code  showing  these 
animations  to  any  interested  reader. 

Summary* 

We  investigate  the  family  of  polynomials  p{z)  —  zn  +  zk  —  1  and  characterize  the  number 
and  location  of  unimodular  roots  in  terms  of  conditions  on  n  and  k*  The  main  result  shows  that 
unimodular  roots  occur  if  and  only  if  n/g  +  k/g  is  divisible  by  6,  where  g  =  gcd(n,  k).  In 
this  case,  there  are  exactly  g  pairs  of  conjugate  unimodular  roots* 
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